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Abstract 

We  study  the  solutions  of  difference  equations  which 
approximate  the  single  nonlinear  conservation  law 


subject  to 


|^+|jf{u)=0  (1) 


u(x,0)  =  u  (x)  ,     u,f  e  R 


Approximations  to  (l)  can  be  calculated  by  explicit  difference 
schemes  due  to  Lax  and  Wendroff  [10,11,15].   We  establish  the 
existence  of  traveling  waves  of  these  difference  schemes  which 

are  approximations  to  shock  waves  of  (l).   That  is,  there  exist 

th 
lattice  functions  for  which  the  n    iterate  of  the  lattice 

function  is  equivalent  to  translation  by  a  distance  (snAt).   The 

speed  s  is  the  same  as  the  speed  of  the  analogous  shock  solution 

of  the  differential  equation. 


V 


1.  Introduction 

¥e  are  interested  in  the  behavior  of  approximations  to  shocks 
which  are  solutions  of 

|^+Jjf(u)=0  (1.1) 

subject  to 

u(x,0)  =  u  (x)  ,  u,  f  e  R 

First,  we  give  a  twenty-five  cent  tour  of  the  theory  of  the  partial 
differential  equation. 

The  nonlinearity  of  (l.l)  introduces  phenomena  which  do  not 
occur  in  the  solutions  of  linear  equations.   In  general,  a  solution 
u(x, t)  will  develop  Jump  discontinuities,  which  are  called  shocks, 
for  initial  conditions  u  (x)  which  are  infinitely  dif ferentiable. 
It  is  necessary  to  extend  the  class  of  solutions  of  (l.l)  to  admit 
these  discontinuous  solutions.   An  integrable  function  u  is  said  to 
be  a  weak  solution  of  (lol)  if  it  satisfies 

00  CO 

//  [^  u(x,t)+^  f(u)]dxdt  +  /   ?//(x,  0)u^(x)dx         (1.2) 

-OO  -CD 

for  all  functions  ?//  which  are  continuously  dif  ferentiable  and  of 
compact  support.   Integration  by  parts  shows  that  a  solution  u(x, t) 
of  (1.2)  satisfies  (l.l)  in  those  regions  where  it  has  one  continuous 
derivative  with  respect  to  each  of  x  and  t.   Suppose  that  u(x, t)  has 
a  jump  discontinuity  which  is  described  by  a  smooth  curve  x(t). 
Moreover,  let  u  (t)  and  u  (t)   be  the  limits  from  left  and  the  right 
of  u(x, t)  at  the  discontinuity.   By  choosing  a  sequence  i   of 


nonnegative  functions  with  supports  which  diminish  to  a  point  on 
the  discontinuity,  the  definition  of  a  weak  solution  implies  the 
Rankine-Hugoniot  relation 

s(u^(t)  -u^(t))  =  f(u^)  -f(up  (1.3) 

where  s  =  dx/dt  is  the  speed  of  propagation  of  the  shock. 

In  general,  there  is  not  a  unique  function  which  satisfies 
(1.2)  for  a  given  initial  condition  u  (x).   The  ambiguity  results 
from  spurious  non-physical  shocks  being  present.   We  give  a  condi- 
tion which  restricts  the  class  of  discontinuities  which  are 
admissible.   We  assume  that  u  <  u  ,  then  we  require  that 

f(u)  <  (u^-u)f(u^)+  (u-u^)f(up  (1.4) 

for  u  in  the  interval  (u  ,  u  ).   Geometrically,  f(u)  lies  below  the 
chord  in  the  (u, f(u))  plane  connecting  the  points  (u  ,f(u  ))  and 
(u  ,f(u  )).   If  u   >  u  the  relation  is 

f(u)  >  (u^-u)f(u^)+  (u-upf(u^)  (1.5) 

for  u  e    (u.,u  ).   Equations  (1.4)  and  (I.5)  are  entropy  conditions 
which  we  abbreviate  to  Condition  E.   They  are  due  to  Oleinik  [23]. 

Hopf,  Lax,  and  Oleinik  [8,9,10,11,20]  have  studied  the  partial 
differential  equation  (l.l)  in  detail.   Recently  Quinn  [22]  has 
established  that  (l.l)  is  well  posed  in  the  t'-.    norm  when  Condition  E 
is  used  to  select  the  physically  relevant  shocks,   Glimm  has 
established  the  existence  of  weak  solutions  for  all  time  when  u  and 
f(.)  are  vectors  in  R  and  u  (x)  has  small  total  variation.   The 


asymptotic  properties  of  these  solutions  are  studied  in  [5] • 
Nlshida  [1?]  has  removed  the  condition  that  u  (x)  be  close  to  a 
constant  for  the  equations  of  isothermal  motion  of  a  gas  in 
Lagrangean  coordinates.   The  twenty-five  cent  tour  of  the  theory  of 
the  partial  differential  equation  is  complete.   The  reader  with 
additional  car  fare  is  advised  to  switch  to  the  recent  survey 
article  of  Lax  [24] , 

We  study  conservative  difference  schemes  which  are,  those 
that  can  be  written  in  the  form 


n+1    n  At 

u  .     =  u  .  --r— 

J       J   Ax 


rg(u.,,,u.,,  -,,....u.  ,  ,-,)-g(u.,,  -,,.o.,U.  ,  )]    (l.Dj 


where  u.  =  u(nAt,jAx)  and  g (•,*,...,•)  is  a  function  of  2k  variables. 
J 

The  equation  (1.6)  will  be  a  consistent  approximation  to  (l.l)  if 
g(u, u, .o.,u)  =  f(u),   A  prototype  of  an  equation  of  the  form  (l.l) 
which  we  previously  discussed  is 

|H.|,f(u)=0.  (1.7) 

Lax  [10]  proposed  the  scheme 

n      n 
,,1^+1    ^J+1  _L   J-1    At  f  „/  n   N   „/  n   >,-,  f-t    q\ 

^  ■  =    — T^ +  —7, -  rTT —  rf(U.,-,)-f(u.-,j]  (loOJ 

J       2      2     2Ax  *■  ^  J+1      J-1 

to  approximate  (1.7)0   This  equation  can  be  written  in  the  form  (1.6) 
taking 

2   2 
,      ,  \    b-a  Ax  ,  a  +b  f-i    r^\ 

g(a,b)=-^.^  +  — ^ —  .  (1.9) 


The  right-hand  side  of  (lo8)  is  a  monotone  increasing  function 

of  each  of  u^,-,  and  u]   -,    so  long  as  l-r-—  u.l  <  1  for  all  j.   The 
J+1      J-1  'AX   J '  ^ 

bound  on  At/Ax  is  the  Courant-Friedrichs-Lewy  stability  condition 

for  the  linearized  difference  scheme. 

In  general,  we  require  that  the  right-hand  side  of  (lo6)  be  a 

monotone  increasing  function  of  each  of  the  variables  u.,,,...,u.  ,  .V 

J+k      j-k 

We  give  the  right-hand  side  of  (1.6)  a  name  rewriting  the  equation  as 

n+1   _,  /  n    n  n   \  /  -,    -,  ^\ 

Note  that  G(u,  u,  ...,u)  =  u. 

We  construct  traveling  wave  solutions  of  (1.6)  with  limits  u 
and  u  at  plus  and  minus  infinity  respectively.   We  require  that  f(«) 
satisfy  Condition  E  across  the  interval  (u  , u  ).   In  the  following, 

-^    Jo 

we  assume  u  <  u„.   The  case  u   >  u„  can  be  reduced  to  this  case  by 
r    £  r    ^ 

symmetry.   A  traveling  wave  u  which  moves  with  speed  s  satisfies 
the  difference  equation 

^x-s(At/Ax)  =  ^^^x+k^'^x+k-l'-'-'^x-k^  "  ^^-^^^ 

The  speed  s  is  related  to  u   and  u.  by  the  Rankine-Hugoniot  relation. 


v-^  That  (1.9)  does  not  depend  explicitly  on  u .  is  a  technical  problem 

J 

that  can  be  removed  by  iterating  the  equation  once.   The  resulting 

relation  is  of  the  form  (1.6)  and  defines  u]        as  a  function  of  u^  „, 

n       n  J  ~ 

u.,  and  u.,„.   Leibniz's  rule  for  calculating  the  partial  derivatives 

of  a  composite  function  implies  that  the  iterated  relation  is  a 

n     n       n 
monotone  increasing  function  of  each  of  u .  „,  u.,  and  u.,„. 


Condition  E  is  the  weakest  condition  that  one  can  expect  to 
imply  the  existence  of  traveling  waves  for  all  values  of  At/Ax. 
Consider  the  closely  related  problem  of  constructing  traveling  waves 
of 

when  e  >  0.   Traveling  waves  exist  and  are  continuous,  Foy  [4] .   To 
obtain  the  relation  they  satisfy  write  u(x, t)  =  v(x-st)  and  substitute 
to  obtain 

-sv'  +  [f(v)]'  =  ev"  . 
Integrate  this  relation  once, 

V'  =  ^^'  ^^^^  +  constant  .  (1.12) 

£ 

The  boundary  conditions  u    =  u  and  u    =  u  determine  the 

I  (JO         X  ~"  LU         JO 

constant  as  that  value  at  which  u^  and  u  are  critical  points 

of  (1.12),  i.eo,  the  right  side  vanishes  when  v  equals  u^  or  u  » 

The  Rankine-Hugoniot  condition  implies  the  right  side  vanishes  at  u^ 

if  and  only  if  it  vanishes  at  u  .   Note  that  no  orbit  running  from 

u   toward  u„  can  get  past  the  first  zero  of  the  right  side  of  (l.ll), 
r         P, 

that  is,  past  a  point  where  Condition  E  is  violated. 

The  traveling  waves  we  construct  are  fixed  points  of  an  opera- 
tor T  defined  by 


(TV)^  =  G^^X+k-n'^X+k-l-Tl'---^X-k-Tl 


) 


where  r)  =  -sAt/Ax.   In  a  subsequent  paper,  we  shall  show  that  these 
fixed  points  are  stable  in  the  following  strong  sense.   If  u^  is  a 
function  defined  on  the  integers  which  takes  on  values  for  which 
G{' ,',...,'  )    is  a  monotone  increasing  function  of  each  of  its  argu- 
ments, then  T%  converges  in  L-j^  to  a  traveling  wave. 


2.  Existence  of  Traveling  Waves  When  ri  is  Rational 

We  fix  a  choice  of  u  and  u.  across  which  f(*)  satisfies 

Condition  E.   The  speed  s  is  determined  from  the  Rankine-Hugoniot 

relation  on  u  ,  u  ,  and  f(*)-   The  parameter  r]   is  then  a  function  of 

6   =   At/Ax.   We  construct  in  this  paper  a  traveling  wave  for  each 

choice  of  6    in  the  interval  (0,0    )  where  6  is  the  largest  value 

^    max         max 

of  6    for  which  each  of  the  partial  derivatives  -^ — (u,,u,  -,,  ...,u_,  ), 
k  >  j  >  -k,  is  nonnegative  when  all  u.,  k  _>  j  _>  -k  lie  in  the  inter- 
val  [u^,u^]  . 

The  construction  of  traveling  waves  is  done  in  two  parts. 
First  we  construct  waves  when  6    is  chosen  so  that  r)  =  s0  is  rational. 
Traveling  waves  which  satisfy  (1.12)  when  s0  is  irrational  will  be 
constructed  as  the  limits  of  traveling  waves  which  satisfy  (1.12) 
with  values  of  s0.  rational  which  converge  to  s6 , 

In  the  remainder  of  this  section,  rj  is  rational.   We  seek 
solutions  of 

^X+Tl  =  ^(^X+k'^X+k-l---^X-k)  ^2.1) 

such  that  u,    =  u  and  u    =  u..   We  iterate  (2.1)  once  and  obtain 
+  oo    r      _  CO     i  ^         ' 
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the   equation 


^x+2ti   "  ^         (^x+k'^^x+k'-l^  •••'^x-k' ^    • 


G^^  is  a  monotone  increasing  function  of  each  of  its  arguments  as 
G  is.  We  may  iterate  (2.1)  a  finite  number  of  times  so  that  N  is 
an  integer,  maintaining  the  form  of  the  equation  after  each  itera- 
tion. Hence,  without  loss  we  consider  (2.1)  when  ri  is  an  integer. 
The  Courant-Friedrichs-Lewy  condition  requires  that  |ti|  <  k.  We  fix 
u  ,  u„,  s  and  f(0  and  consider  functions  u^  defined  on  the  integers. 
We  seek  traveling  waves  as  solutions  of  (2.1)  with  limits  at  plus 

and  minus  infinity,  ^^  qq  =  ^-^   ^i^^  ^_  od  ^  ^H' 
We  prove  the  following 

Theorem  1:  {t\   rational)   Let  f(')  satisfy  Condition  E  across 

fu  ,u  )  where  u  and  u„  determine  a  shock  solution  of  (l.l)  which 
^    r      S,'  r      i 

travels  with  speed  s.   Suppose  that  Q    =   At/Ax  is  chosen  so  that 
sAt/Ax  is  rational  and  is  smaller  in  magnitude  than  k.   Moreover,  we 
suppose  that  the  right-hand  side  of  (2.2)  is  a  monotone  increasing 
function  of  each  of  u  _^,,-u.   +v_i' •  •  •  ^  ^x-k  ^^^"^  ^^^  arguments  take  on 
value  in  the  interval  [u  jU  ].   Then  for  each  number  u^  e  {\i^,u^), 
there  is  a  unique  function  defined  on  the  integers  which  takes  on  the 

value  u  at  X  =  0  and  which  satisfies 
o 

^X+Tl  =  ^X  -  1^  fg(^X+k'^X+k-r'--^X-k+l^ 

(2.2) 
-  S^^x+k-r^x+k-2"--'^x-k)^  • 

The  solution  u   is  a  monotone  function  of  x  and  depends  continuously 
at  each  value  of  x  on  u  . 
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The  monotonicity  of  the  right  side  of  (2o2),  which  we  have 
named  G( •,•,...,*)  previously,  is  central  to  the  proofs  which  follow 
We  collect  the  consequences  of  this  assumption  that  we  will  need 
later.   We  name  the  variables  on  which  G  depends  writing  it  as 
G(u,,u,  -,,...,  u  ,  ).   The  partial  derivatives  with  respect  to  each  of 


u,,,u,  .,  ...,u  ,  are  positive.   Hence,  writing  g  =  g(^j^^^i^_]_^  •  •  • ' w_1^4-]_) 


and  4^  =  gi  >  it  follows  that 
dw^   ^k' 

5  .  -  e(g.  -g.,  J  >  0   for   k  >  j  >  -k 

J  J       J   -^ 

where  5  .  =  1  if  j  =  0,  and  0  otherwise.   Moreover,  gi,^-|  =  S  k  "  *-** 
From  g,  ,  T  =0,  it  follows  that 

-eg^  >   0 

and  by  induction,  -6g.    >  -6g.,n,    3    >   0   implies 

-eg.  >  0   for   k  _>  j  >  0  .  (gl) 

At  j  =  0,  we  obtain  (l  -0g  )  >  -6g,      >   0.   moreover,  by  induction  it 
follows  that 

(1-  eg.)  >  0   for   -k  <  j  <  0  .  (g2) 

J 

Additional  inequalities  can  be  obtained  by  considering  first  that 
0g  ,  >  0  and  then  proceeding  by  induction  for  j  >  -k.   We  need  only 
the  relations  (gl)  and  (g2). 

¥e  give  now  an  estimate  for  the  difference 

^(^x+k'^x+k-l-'-^x-k)  -^^^x+k^^x+k-l---^x-k)-   ^^  ^^^  ^^^  "'^^'^ 
value  theorem  and  write  the  difference  as 
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/  5i:|^(P<Mi-e)v;)K,j-v^,j)dp 

Q   J — i^     J 

Where  u*  =  (^x+k' ^x+k-1' ' *  * ' ^x-k^ *   ^^®  partial  derivatives  are 
positive.   Hence  the  absolute  value  of  the  difference  is  the  above 
expression  but  with  ( u^_j_  .  -  v^_^  )  replaced  by  |  u^_j_ .  -  v^_^  .  |  .   Using 
the  mean  value  theorem  again,  we  get 

^(^x+k'^x+k-r  •  • -^x-k)  -  ^^^x+k'^x+k-r  •  • -^x-k) 

(p) 

<  |u^-v^|-0Vg(7^)-|u^-v^l+eVg(7^_l).|u^_l-v^_ll 

where  -y   is  a  point  on  the  chord  joining  ^^  ^  ^^x+k' ^x+k-1' *  *  * '  ^x-k-1^ 
^^^  ^^x+k'^x+k-1' *  *  "^x-k+l^*   Inequality  (P)  is  an  equality  if  and 
only  if  each  of  u^^^  -  v^^^,  ^x+k-1  "  ^x+k-1' '  * '  ^"^^  ^x-k  "  ^x-k  ^^^  ^""^ 
same  sign. 

Prior  to  proving  that  a  traveling  wave  exists,  we  prove  that, 
assuming  existence,  the  solution  must  be  a  monotone  function  of  x 
and  that  a  solution  depends  uniquely  on  the  value  which  it  assumes 
at  a  single  point.   Subsequent  to  these  proofs  we  construct 
traveling  waves. 

Lemma  2.1.   Let  u  and  v  be  two  traveling  waves,  such  that 

—  J^  x 

u  ,v   e  [u  ,uj  for  all  x,  then  u  and  v   agree  at  one  point  of  the 
X  X     r   iJ  XX 

lattice  if  and  only  if  they  agree  at  every  point.   If  u^  >  v^  at  a 
single  point,  then  u^  >  v   for  all  x. 

Proof:   By  assumption 


u  ,   —  Giu  .-ijU  ,,  -i,o..,u  , 
x+T]     ^  x+k'  x+k-1'    '    x-k' 


and 


then  by  inequality  (P) 

Sum  this  inequality  from  x  =  -M  to  +M,  and  eliminate  the  common  terms. 


,  /^   /^ 


As  M  — ►  +  00,  each  of  the  terms  goes  to  zero.   We  obtain 

0  <  0  . 

This  implies  that  Inequality  (P)  was  an  equality  for  each  value  of  x. 
Hence  v^^^  -  u^^^,  u^^j^_^  -  v^^^_^,  . . .  and  v^_^  -  u^_j^  are  of  the  same 
sign  for  each  value  of  x.   The  sign  may  depend  on  x.   The  difference 
V  -u  cannot  go  from  positive  to  negative  without  being  zero.   If 
u^  =  V  for  a  single  value  of  x,  then 

X      X 

1 

0  =  u-v  =  r  |S  (eu*   +(l-0)v*   ).(u*-v*)d0 

X    X    J    dj  ^   X-TJ   ^        X-T]      X    x' 

0 

but,  as  ^G/^j  >  0,  it  must  be  that  u,  =  v,  for  t  e    [x-i^-k,  x-ri+k]  . 
Proceeding  by  induction,  it  follows  that  u  =   v   if  u  =  v   at  a 

XX        XX 

single  point.   Thus  if  u^  >  v^  at  a  single  point,  u^  >  v   for  all  x. 

XX  A       A 

This  completes  the  proof. 

Consider  that  u  is  a  traveling  wave,  then  u  -,  is  also  a 
traveling  wave.   Hence,  u^  -u^,-,  is  of  the  same  sign  for  all  x. 

X     X~r  X 

Therefore  we  have 
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Lemma  2.2.   Let  u  be  a  traveling  wave  such  that  u  e  [u  ,u.], 
X  X     r   i; 

then  u   is  a  monotone  function  of  x  and  lies  in  the  interval  (u  ,u„). 
X  ^    T      i' 

¥e  now  proceed  to  the  construction  of  traveling  waves.   ¥e 
study  the  fixed  points  of  the  operator  T   defined  on  functions  whose 
domain  is  the  integers  greater  than  -(k+r)).   For  x+r|  >  0  we  define 

^V^x+Ti  =  ^^^x+k'^x+k-l---^x-k)  (2.3) 

T  is  the  identity  on  the  integers  x  such  that  -k-  t]  <  x  j<  0.  We 
define  a  fixed  point  of  T  to  be  a  function  u  defined  on  1  ,  the 
integers  greater  than  -k-T],  which  satisfies 

V   =  (T  v) 
X    ^  r  'x 


and  has  the  limit  u,    =  u  . 

+  00    r 

A  one-parameter  feunily  of  fixed  points  of  T   can  be  extended 
to  a  traveling  wave.   We  identify  this  family  as  those  fixed  points 
of  T  which  can  be  "patched"  to  a  fixed  point  of  T  .   T   is  analogous 

-'■  Jo  X/ 

to  T   and  is  defined  on  functions  with  domain,  1~ ,    the  integers  less 


than  k-T]  by 


T.w)  ,   =  G(w  ,,,w  ,,  -,,...,w  ,  ) 
■  i    'x+T)     ^  x+k  x+k-1    '    x-k' 


when  x+T]  <  0  and  (T  w)  =  w  when  0  <  x  <  k-ri.   We  require  a  fixed 

i/   X      X  — 

point  of  T.  to  have  the  limit  w    =  u„.   We  define  the  initial  value 
i  -  CO     i 

of  a  function  with  domain  equal  to  I   to  be  the  values  which  it 

assumes  on  the  interval  -k  -  ri  <  x  _<  0.   Here  and  in  the  following, 

X  assumes  values  in  the  integers.   The  objective  of  the  following 
lemmas  is  to  establish 
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Lemma  2.3'   Let  ^   be  an  initial  value  on  I   such  that 
^r  —  ^x  —  ^9    ^°^   -k  -  T)  <  X  <  0  and  u  _<  |   <  u  ,  for  x  =  0,  then 
there  exists  a  unique  fixed  point  u  of  T  with  ^   as  initial  value. 
For  a  fixed  value  of  x,  u  is  a  continuous  function  of  |.   Moreover, 

X 

u  <  u  <  max  P      for  x  >  0. 
r  —  X       ^x 

An  analogous  result  is  true  for  T,o   The  initial  value  of  a 
function  defined  in  I",  is  its  value  on  0  j<  x  <  k-r^.   We  state 

Lemma  2.4.   Let  C   t)e  an  initial  value  on  I~  such  that 

u   <  C   ^  u.  for  0  <  X  <  k-ri  with  u   <  C   ^  u.,  at  x  =  0,  then  there 

exists  a  unique  fixed  point  w^  of  T.  with  C^  a-s  initial  value.   For 

^  &  X 

a  fixed  value  of  x,  w   is  a  continuous  function  of  C*   Moreover 

X 

min  C^  <  w^  <  u^. 

With  these  two  lemmas  we  can  construct  traveling  waves  when  r\ 

is  rational.   We  prove  Theorem  1  now  assuming  Lemmas  2.3  and  2.4. 

Fix  u^  e  (u^,u^)  and  let  |^  satisfy  u^  <  |^  <  u^,  -k  -  ti  <  x  <  0, 

with  I   at  X  =  0  equal  to  u  .  Then  there  exists  a  fixed  point  of  T  , 

u^d)  such  that  u   <  u  <  u.  for  all  x.   The  function  u  restricted 
X  r  —  X     j>  X 

to  the  interval  [0,  k-ri)  is  an  initial  value  of  a  fixed  point  w^  of  T.. 

Moreover,  ^y-   "^  "^x  —  ^P    ^'~^^   ^^'^   "^  ^"^"^  -^^^  -k  -  r)  <  x  <  0  in  particular. 
The  fixed  point  w  determined  in  this  manner  is  a  continuous  function 

X 

of  ^  .   Moreover,  for  all  ^,  u  <  2   <  u,  for  -k  -  ri  <  x  <  0,  it 
X  r  —  X  —  a, 

follows  that  u  _<  w  (^)  <   u   for  -k-  t]  <  x  <  0.   Hence  the  Brouwer 
fixed  point  theorem  implies  there  is  a  fixed  point  |  .   Define  a 
function  v   on  I   as  u(|*)  and  on  I~  as  w  {i* ) .      As  vi  (^*)  =  ^*,  the 

XX  XX 

definition  is  the  same  on  I  n I~ .   Notice  that  v   is  a  traveling  wave 

X 

of  the  difference  scheme  and  takes  on  the  value  u  at  x  =  0,  which 
was  chosen  arbitrarily.   We  gather  this  into 
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Lemma  2.3 ■   Theorem  1  is  a  consequence  of  Lemmas  2.3  and 
2.4. 

Proof:   All  that  remains  is  to  show  that  a  traveling  wave 
depends  continuously  on  the  value  it  assumes  at  any  point. 

Let  V  be  the  traveling  wave  which  takes  on  the  value  u.  at 
X  =  0  and  suppose  that  u.  -+  u  .   We  may  suppose  that  u   <  u.  <  u  for 
all  i  and  that  u   and  u~  e  (u  >u  ).   Then  by  Lemma  2.1, 

-     i    + 

u   <  V   <  u 

XXX 

-       +  -       + 

where  u  and  u  are  traveling  waves  which  u   and  u  determine 

X         X 

respectively.   Extract  a  subsequence,  if  necessary,  so  that  v"""  — ^  v  . 
Then 

^x+T  =  ^^^x+k'^x+k-r-'-^x-k) 

and  moreover,  u=u7    "^v,    <u,    =u.   Similarly,  v    =  u„. 
r    +00—  +00—  +00    r  "^-oo    I 

This  finishes  the  proof. 

We  now  turn  to  the  proof  of  these  Lemmas.   We  assume  that  all 
functions  defined  on  I  take  on  their  values  in  the  interval  Tu  ,uj. 

First,  we  show  that  a  fixed  point  of  T   depends  uniquely  on  its 
initial  value.   This  we  do  directly  from  inequalities  (gl)  and  (g2). 
To  establish  the  existence  of  fixed  points,  we  consider  the  set  (X  of 
a  e  [u  , u  ]  such  that  every  initial  value  for  which 

^v  "^  ^   ^  ^   when   X  <  0 


and 


u   <  I   <  a    for   X  =  0 

X  —^         X 
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determines  a  fixed  point  of  T  .   We  show  that  Ci   is  both  open  and 
closed  in  the  Interval  [u  ,  uj.   Using  a  continuity  argument. 
Lemma  2.5  follows  from  the  existence  of  fixed  points  with  initial 
values  close  to  u  .   Lemma  2.4  then  follows  by  symmetry. 

Lemma  2.6.   (Unique  Dependence  on  Initial  Values)   Let  v   and 

w  be  fixed  points  of  T  ,  which  have  the  same  initial  value,  then  v 
X  -^  r'  '       X 

is  equal  to  w  for  all  x. 

Proof:   Difference  the  relations  that  v   and  w   satisfy,  apply 
inequality  (P),  and  sum  from  x+r)  >  0  to  x+r]   =   M.   Eliminate  the 
common  terms  from  the  sums  on  each  side  of  the  inequality  and  let 
M  -^  +  OD.   The  result  is 

x<k-'q 

V  -  w  I  -  evgf-Y  )•  |v  -  W  I  <  0 

X    x'      '=''''o'''o    o'  — 

x>( 

where  -v   is  a  point  in  R   on  the  chord  between  v^  and  w^.   Using  the 
<  o  o      o 

assumption  that  v^  =  w   for  x  <  0,  this  becomes 

X  X  — ■ 

X<k-T] 


where 


^      °^x-l^x-^xl    1  0 
x>0 


^      -   ^     c^U     ^^  ^-Yo)  ^°^  k-Tl     >    X     >     -Tl 


a. 


X 


x+r)+l 


X+Tj+l 


That  5^  >  0  follows  from  (gl)  and  ( g2 )  and  hence  v   =  w   for 

X  XX 
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0  <  X  <  k-T).   Translation  to  the  left  on  I  maps  fixed  points  into 
fixed  points.   Thus  v^  =  w^  for  x  >  0  by  induction. 

We  now  give  an  estimate  on  the  values  of  a  fixed  point  in 
terros  of  the  initial  value  .   This  is  a  maximum  principle  and  does  not 
require  that  the  lattice  function  v  has  a  limit  at  x  =  +  od. 

Lemma  2.7.   (Maxim\m  Principle)   Let  v   satisfy 

^X+Tl  =    G^^X+k'^X+k-l"-"^X-k)     ^°^    ^^^  "    ° 

V  e  fu  ,u  )  for  all  x,  then  v   assumes  its  maximum  (minimum)  at  a 
X    ^  r   ^  X 

value  X   >  0,  if  and  only  if  v   is  a  constant.   No  condition  is 

O  X 

imposed  on  the  behavior  of  v   as  x  — ►  +00. 
Proof:   For  x   >  0, 

^X   =  ^^^X  +k-n'%  +k-Tl-l---^X  -k-Tl)   •  ^2.4) 

o       o    '   o    '         o 

Moreover,  assume  without  loss  that  v^  >  v   for  all  x.   G  is  a  strictly 

x^-     X 

monotone    increasing   function   of  each  of   the    arguments.      Hence 

V         <   G(v      , V      , . . . , V      )    =  V 
X      —       ^    X         X  '     X    '  X 

o  o  O  DO 

and  the  inequality  is  strict  unless  v   _^+-j  =  ^   for  -k  j^  j  _^  k. 

o  ^  "^     o 
But  (2.4)  is  a  (nonlinear)  recursion.   If  a  solution  is  constant  for 

any  2k  consecutive  values,  the  solution  is  constant  for  all  x.   This 

completes  the  proof. 

We  prove  a  weak  result  on  the  existence  of  fixed  points.   The 

result  implies  that  the  same  inequality  on  each  component  on  two 

initial  values  persists  to  the  traveling  waves. 
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Lemma  2.8.   (Partial  Ordering  of  Fixed  Points)   Let  u  =  ^   for 

-k-T]  <  X  <  0  be  the  initial  value  of  a  fixed  point  u  o   For  each 

initial  value  t    ,    u      <   t      <   t    ,    there  is  a  fixed  point  w   of  T  which 
^x'   r  —  ^x  —  ^x  ^  X     r 

takes  on  the  initial  value  ^      and  w  <  v   for  all  x. 

Proof:   By  the  maximum  principle,  u   >  u  for  x  >  0,  unless 

t      =  u  .   Therefore,  the  lattice  function 
^x    r 


o 
w 

X 


r  C     if   X  <  0 
^x  — 


u    otherwise 
r 


is  bounded  from  above  by  u  .   Define  w   to  be  the  n   iterate  of  w 

'^    X  X  X 

under  T  .   As  G  is  a  monotone  increasing  function  of  each  of  its 
r 

arguments 

^x+T  =  ^^^S+k^^x+k-r-'-'^x-k^  -  ^r 

for  x+T)  >  0.   Moreover,  by  induction 

_^  n  _^  n+1  ^ 
u   <  w  <  w    <  u   . 
r  —  X  —  X   —  X 

Thus  the  limit  w  as  n  — ►  co  exists  which  we  write  as  w  and  satisfies 

X  X 

u   <  w  <  u 
r  —  X  —  X 


Hence  w,    =  u  ,  and  the  proof  is  complete. 
+  00    r'         ^  ^ 

The  continuous  dependence  of  a  fixed  point  upon  its  initial 
value  now  follows. 

Lemma  2.9.   (Continuous  Dependence  on  Initial  Values)   Let  ^ 
be  a  sequence  of  initial  values  such  that 
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r  —  ^x  —  ^x 

where  C   is  an  initial  value  that  determines  a  fixed  point  of  T  , 
then  if  l"^  -►  |   for  each  value  of  x,  it  follows  that  v  (^"'')  ->■  v  (O 

XX  XX 

pointwise,  where  v  i^'^)    is  the  fixed  point  of  T  that  ^   determines. 

X  -^ 

Proof:   Because  of  the  monotonicity  of  fixed  points  of  T^,  it 
follows  that 

Passing  to  a  subsequence  if  necessary,  it  follows  that  v  (^  )  — ►  w  . 

X  X 

But 

from  which  it  follows  that  w    =  u  o   The  uniqueness  result  implies 

We  define  (f(a)  to  be  a  set  of  initial  values,  ^  ,  such  that 

u  <  ^   <  a   for    -k-Ti  <  x  <  0 
r  —  ^x  —  ' 


and 


u     <   i      <  a   when    x  =  0 
r  —  ^x 


We  are  interested  in  the  supremum  of  the  set  <^  of  a  e  [u  ,u  ]  such 

that  every  initial  value  in  u-(a)  admits  a  fixed  point  of  T  .  Call 

this  point  a*.   Then  Cl  ^    [u^,  a*)  or  [u^,  a*]  by  Lemma  2.9.   We  first 
show  that  (^   is  closed. 

Lemma  2.10.   (  Cl    is  closed)   Let  a*  and  cl   be  as  above  then 
every  initial  vaj.ue  in  /(a*)  determines  a  fixed  point  of  T^. 


17 


Proo 


f:   Let  a.  t  a*  where  a.  e  C'L .      Suppose  P      is  an  initial 


value  in  /(a*).   Hence  u  <  |   <  a*  and  we  may  suppose  that  ^   =  a^ 
for  some  x  e  (-k-rj,  O).   We  proceed  as  before  setting 


o 
u 

X 


^'  P  if        X   <   0 

^x  — 


u  if        X   >   0 

1^      r 


n        n^n  o 


Then   u      =   Tutu   }    which   satisfies 
X  r  X  '     x^ 

^x+r|   =   ^^^x+k'^'x+k-l'-'-^^x-k^    • 

Moreover  as  ?   <  a  and  ^   =  a*  for  some  x  the  maximum  principle 

implies  that  u   <  a*  for  all  x  >  0  and  thus  u   restricted  to  (0,k+Ti] 

is  an  initial  value  ^  which  lies  in  (^(a  )  for  i  sufficiently  large. 

Hence,  there  exists  a  fixed  point  of  T^,  w  ,  with  initial  value  C« 

r   X 

'^  o  n 

Set  w  =  w   /i  ,  \»   Then  w   >  u  for  x  >  0  and  as  w  bounds  u  on 

X      X—  y  K.~rT1  j  X       X  XX 

the  interval  (O,  k+ri]  it  follows  that  u^  <  w   for  x  >  0.   Hence 

X  -~   X 

u   <  u  <  w   and  u,    =  u  .   This  completes  the  proof, 
r  —  X  —  X      +00    r  i-  jr^ 

We  now  turn  to  characterizing  a*. 

Lemma  2.11.   Let  a  >  u  .   If,  there  exists  a  fixed  point  of  T 
with  initial  value  |   =  a,  then  {^i   contains  an  interval  about  a. 

Proof:   We  assume  that  there  is  a  fixed  point  u  of  T  with 
initial  value  ^^  =  a.   As  u   is  not  a  constant,  it  follows  from  the 

X  X 

maximum  principle  that  u^  <  a  for  x  >  0.   The  restriction  of  u  to 

X  X 

(0,k+ri]  is  an  initial  value  ^  which  lies  in  (f(a-e)  for  some  e 
positive  and  sufficiently  small.   The  dimension  of  (5(a-e)  is  k+rjo 
Consider  the  equation 
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^x+Ti  =  ^(^x+k'^x+k-r---^x-k 


) 


as  a  nonlinear  recursion  of  order  2k.   The  process  of  calculating 


the  value  of  v  ,  as  a  function  of  "^j^+j^^  ^x+k-1' "  *  * '"^x-k+1  ■^^  ^ 

2k    2k 
function  Z.'  in  the  usual  way  from  R   ->  R  ,  which  is  well  defined 

locally  so  long  as  ^—    (v)  7^  0.   This  makes  -^  an  open  mapping.   The 

du^ 

continuous  dependence  of  a  fixed  point  on  its  initial  value  implies 
that  the  set  of  vectors  in  R   near  {y-^,v^,  . . ,  ,v^^)   which  generates 
fixed  points  of  T   is  a  manifold  of  dimension  k+r).   Call  this  set  ')?!!. 
Then  X'     '^{■]ri)    is  of  dimension  k+ri,  because  'Z'   is  an  open  mapping. 
Fixed  points  of  T  depend  uniquely  on  their  first  k+rj  components. 
The  projection  of  ■^'^'^'^ {'))?)    onto  its  first  K+t]   components  must  be  of 
dimension  k+rj.   This  completes  the  proof. 

If  we  suppose  that  a*  is  greater  than  u^,  this  lemma  enables  us 
to  determine  a*.      Select  a  sequence  a.  f  a'^  then  the  fixed  points  u^ 
such  that  u^  =   a}   for  -k-r)  <  x  <  0  satisfy 

for  all  X  and  i.   Hence  u?;  f  w^  at  each  value  of  x.   If  w  is  not  a 
constant  then  w^  <  a*  for  x  >  0  and  as  before  w^  is  a  fixed  point 


of  T^. 


Each  of  u"^  is  a  fixed  point  of  T„.   Sum  the  relation 


X 


^i+n  =  ^^^i+k'^i+k-r"--Vk) 


x+r\ 

from  X   >  0  to  X  =  +M,    and   remove   the    terms   in   common  on  both  sides   of 
the   relation.      Let  M-++00.      It   then  follows   that 
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<k-Ti      .         <k- 


y^  u^  =  y^  u"^  -  eg  ( u,^     . ,  11?"     ,,...,  u  ,     )  (2.5) 

X  >  0  x>-'n  '  '  ' 


-0g(u^,  u^,  ...,u^)  +  riu^    . 


Note  g(u,  u,  ...,u)  =  f(u)  and  let  u  -►  w  .   As  w  =  a*,  it  follows 

XX  X 

that 

Ti(a*-u^)    =   0[f(u^)-f(a*)]     . 

Recall  r\   =  -s6   and   u     ^  ^-^  j^  u..      Then   a*    is    the    solution   of 

s(a*-u^)    =   f(a*)-f(v^)    . 

Condition  E  implies  that  u.  is  the  first  solution  of  this  equation 
greater  than  u  . 

By  the  continuity  method  we  have  proved 

Lemma  2.12.   If  there  exists  an  e  >  0  such  that  each  initial 
value  in  .-'■'(u  +e)  determines  a  fixed  point  of  T  ,  then  every  initial 
value  in  ^{u.)    determines  a  fixed  point  of  T  . 

We  now  turn  to  showing  the  existence  of  fixed  points  of  T  with 
Initial  values  close  to  u  c   We  do  this  in  stages. 

To  begin,  we  strengthen  Condition  E  by  imposing  the  stronger 
shock  condition 

f'(u^)  <  s  <  f'(up  . 

This  relation  with  (^)  replacing  (<)  is  a  result  of  Condition  E.   We 
will  later  weaken  the  strict  Inequality  to  (<)  by  taking  a  limit. 

We  now  establish 

Lemma  2 . I3 .   If  s  >  f ' (u  ),  then  it  is  possible  to  choose  e  >  0 
so  that  each  initial  value  ?   for  which 

X 
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r    _  ^x  —  r 

determines  a  fixed  point  of  T  . 

Proof:   Define  the  set  jP  „  to  be  the  set  of  lattice  functions 
with  domain  I  ,    which  satisfy 


-  u   <  ae  '^ 


'  X   r '  — 

where  a,  p  >  0.   Let  v     e     -f       ,    then  it  follows  from  the  mean  value 
theorem  that 

/„   \    Pfx+Ti)   <  ^     /      \  -Pj+pTi      3(x+j) 
(T  V )  ,  e'^^   '^  =  >     c.(v  ,u  je  ^"^    ^'v  ..e"^^   ^^  , 
^  r  x+Tj  'i_  -.       j^  x'  r'  x+j 

The  c.'s  are  positive  as  a  result  of  the  assumption  of  the  monotonicity 
of  G( •,•,...,•) •   Applying  the  triangle  inequality  under  the  assump- 
tion that  V  £  f-   o  we  obtain 

(T  v)  ^  e^^^+^)   <  ah(a,P) 
where 

h(a,p)  =     sup       rm  ^•(^x'^r)^"^^'"'^'^  )  • 
Cv^:|v^-u^|<a}V  J^^  ^'^ 

An  easy  calculation  yields  that 

h(0,P)  =  (l-eP)eP^  [1-0  XZZ  |g-  (V^r-'-'^r)^"^"^' 

A  Taylor  series  expansion  in  ^  gives 

h(0,P)  =  1-  e(s-  f'(u^))p  +00^)  . 
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We  have  used  f'(u)  =  )    -^^  (u,  ...,u),  which  follows  from 

f(u)  =  g(u,  u,  .  . .,  u) .   Hence,  as  0(s-  f'(u  )  >  0  we  may  choose  a 

P  >  0  so  that  h(0, p)  <  1.   Fix  this  choice  of  p  and  choose  a  >  0,  so 

that  h(a, p)  <  1.   For  this  choice  of  a  and  P,  we  have  shown  that  T 

maps  ^     o  into  itself.  J-      o  is  convex  and  compact  in  the  topology 

of  uniform  convergence.   T  has  a  fixed  point  by  the  Schauder  fixed 

point  theorem.   The  lemma  is  proved. 

Lemma  2.3  is  established  under  the  additional  constraint 
f ' (u  )  <   s.   Lemma  2.4  is  true  by  symmetry.   The  analogous  additional 
constraint  is  s  <  f'(u.).   Theorem  1  is  established  for  every  choice 
of  f (• )  so  that 

f'(u^)  <  s  <  f'(u^)  . 

We  now  relax  the  inequalities.   We  consider  the  traveling 
waves  to  be  fixed  points  of  an  operator  T  which  depends  on  g.   Note 
g(u, u, ...,u)  =  f(u).   We  abuse  the  notation  and  write  T  =  T(f). 

Suppose  that  f  satisfies  Condition  E  and  the  inequalities 

f'(u^)  <   s  <  f(up  . 

Choose  (t  positive  and  smooth  on  [u  ,uJ  so  that  h     vanishes  on  the 
complement  of  (u  ,u  +e).   Set  f  =  f+  ^      and 

2k 
g^{u-^,u^,  ..  .,u^^)    =  g{u-^,\i^,  ...,u^^)    +  >    (t)^(u  )/2k  . 

Furthermore,  choose  cj)   so  that  G  (•,♦,...,•)  is  monotone, 

f^(u^)  <  s  <  f'(u^)  , 
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and  so  that  f  -*  f  on  [u  ,u  ].   Fix  u  e  (u  , u  ).   There  exist 

functions  u^  so  that  u^  =  u  and 
X  o     o 

T(f  )u^  =  u^  . 
^  e  ^  X    X 

Moreover,  there  is  a  fixed  point  w   of  T,(f  )  with  initial  value 
identically  equal  to  u  which  is  independent  of  e  as  T.(f  )  =  T,(f) 
for  lattice  functions  whose  range  lies  above  u  +e .   Hence 

w^  <  uj  <  u^  (2.6) 

for  X  j<  0.   As  u  e  [u  ,u  ]  pass  to  a  subsequence,  if  necessary,  so 
that  for  each  value  of  x,  u^  -^  u^.   Then  T(f )u  =  u  .   As  u^  is 

XX  XXX 

monotone,  u   is  also.   The  limits  u   ,  u    exist;  u_    equals  u 
by  (2.6).   Moreover,  we  obtain 

s(u^   -u J  =  f (u^   )  -  f (u.) 
^+00   ^'  ^+oo^    ^    H 

as  in  the  proof  of  Lemma  2.12.   Hence  u,    =  u  . 
^  +  00    r 

We  now  have 

Lemma  2.l4.   Theorem  1  is  established  for  every  f(*)  which 
satisfies  Condition  E  across  [u  ,uJ  and 

f'(u^)  <  s  <  f'(u^)  . 

We  now  relax  the  last  inequality. 

Lemma  2.15.   If  s  >  f'(u  )  and  there  exists  a  value  u   so  that 
f(«)  satisfies  Condition  E  across  [u  ,  u  ]  then  it  is  possible  to 
choose  e  >  0,  so  that  every  initial  value  in  <5(u  +  e)  determines  a 
fixed  point  of  T  . 
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Proof:   It  is  possible  to  change  f(')  to  f  near  u,  so  that 
fp.  satisfies  Condition  E  across  [u  ^u.]  and  so  that 

f;(up  >  s  . 

Then  there  exists  a  fixed  point  of  T(f  )  which  takes  on  the  value 
(u  +  u  )/2  at  X  =  0.   This  traveling  wave  is  a  fixed  point  of  T  (f) 
on  I  .   Lemma  2.8  completes  the  proof  of  this  lemma. 

Moreover,  Lemma  2.3  is  proved  in  complete  generality.   Lemma 
2.4  is  true  by  symmetry.   Theorem  1  is  established  for  t]  rational. 
In  the  next  section  we  show  these  traveling  waves  are  stable  as  a 
sequence  tj  .  of  rational  numbers  converges  to  an  irrational  number. 


3 .  T)  Irrational 

A  good  difference  scheme  for  calculating  approximations  to 

discontinuous  solutions  of  a  partial  differential  equation  produces 

approximations  which  change  from  values  near  u.  to  values  near  u  in 

a  few  mesh  widths.   Experience  obtained  from  numerical  experiments 

indicates  the  behavior  of  the  transition  between  u„  and  u  for 

^      r 

traveling  waves  is  typical  of  the  behavior  of  approximations  to  shocks 

generated  by  the  finite  difference  scheme  for  an  arbitrary  initial 

condition.   The  results  of  Theorem  1  imply  that  the  traveling  wave 

u  (u  )  is  a  function  of  x  and  u  ,  the  value  of  u   at  x  =  0.   We 
x^  o'  o'  x 

produce  this  function  explicitly  for  the  special  equation  (1.9)  with 
u  =  -1  and  u.  =  +1.  The  speed  s  of  the  traveling  wave  is  zero,  and 
the  equation  for  the  traveling  wave  is 
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^x+l"^  ^x-1    e    ,2  2      s  ,^  -  . 

Equation  (2.5)  specialized  to  this  equation  is 

-    o    (^^4-1  +^v)  +^  •  (3-2) 


x+1 


2     2  ^  x+1   X 

This  is  a  quadratic  relation  which  can  be  solved  explicitly  for  u, 

as  a  function  of  u  .   We  set  u  =  0  and  solve  for  u-, 

X  o  1 


u^(e)  =  +  i^4^  •  (3.3) 

We  choose  the  root  so  that  u-,(e)  <  0.   The  Courant-Friedrichs-Lewy 
condition  for  this  equation  is  0  _<  -p.   As  0  varies  from  0  to  -^, 
u-,(e)  varies  from  0  to 

-  /^+1 

and  is  monotone  in  0o   This  substantiates  the  empirical  observation 
that  the  profile  of  a  traveling  wave  becomes  steeper  as  6   is  increased 
to  the  limit  imposed  by  the  Courant-Friedrichs-Lewy  condition.   We  do 
not  establish  this  result  in  general.   We  establish  a  weaker  result. 
We  show  that  the  traveling  waves  we  have  constructed  converge  to  a 
monotone  profile  as  v\.    runs  through  rational  values  to  an  irrational 
limit.   The  limiting  shape  is  defined  on  the  entire  real  line  and 
satisfies 

^x+n  =  ^(^x+k^^x+k-i---^x-k^  ^^-"^^ 

and  has  limits  u,    =  u  ,  and  u    =  u^. 

+  00     t'  -  go     H, 
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When  T^  is  rational  the  minimal  domain  over  which  O-^)   makes 
sense  is  the  linear  span  over  the  integers  of  1  and  t].      The  distance 
between  adjacent  points  in  the  lattice  is  a  function  of  r|,  which  we 
call  A(t]).   As  ri  converges  to  an  irrational  number  through  rational 
values,  A(ri)  -*   0.   Let  u  (u  ;ri.)  be  the  traveling  wave  which  satis- 
fies (3.4)  with  T]  =  T].  and  takes  on  the  value  u   at  x  =  Oo   We  will 
construct  a  traveling  wave  as  the  limit  of  u^  =  ^x^^o'^i^  ^^  ^i  ~*  ^° 
Fix  r\    irrational  and  suppose  rj^.  rational,  converges  to  i].      Choose  a 
sequence  s.  such  that  ri./s.  =  0  is  constant  and  such  that 

s.(uj  -u^)  =  f(up  -f(u^) 

and  f  satisfies  Condition  E  across  [u^,up.   If  necessary,  f  can  be 
extended  on  [u,,u,+e].   We  define  v^  on  the  entire  real  line  as 
piecewise  constant.   It  is  a  consequence  of  the  partial  ordering  and 
the  maximum  principle  that  a  fixed  point  u^  of  T^  with  r)  =  ri  <   t] 
for  all  i  and  which  satisfies  u^^  " -^o   satisfies 

u,  <uj<u^.  (5.5) 

Extract  a  subsequence,  if  necessary,  so  that  u  -^  u  for  all 

rational  x.   The  limit  u  is  a  monotone  function  and  we  may  assume, 

without  loss  that  u  is  continuous  from  the  left.   Moreover,  an 

integral  fonn  of  (3.4)  and  the  bounded  convergence  theorem  implies 

u  satisfies  (3.4)  for  all  x.   The  inequality  (3.5)  implies  that 

u    =  u  .   As  u   =  u   at  X  =  0,  it  follows  as  in  Lemma  2.12  that 
+  CO     r       X    o 

u     =  u  .  o 

-  00    e. 

We  have  established  the  existence  of  a  traveling  wave  when  r\ 

is  irrational  which  takes  on  the  value  u   at  x  =  0.   We  state  this  as 
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Theorem  2.   (-q  irrational)   Under  the  same  condition  as  Theorem 
1,  except  that  x]   is  irrational,  there  exists  a  traveling  wave  v^ 
which  is  monotone  and  continuous.   Any  other  continuous  traveling 
wave  is  a  translate  of  v  . 

Proof:  To  obtain  the  uniqueness,  systematically  substitute 
integral  signs  for  the  summation  signs  in  the  proof  of  Lemma  2.1. 
We  proceed  to  show  that  the  traveling  wave  constructed  above  is  a 
continuous  function  of  x. 

Suppose  V   and  w  are  solutions  of  (3.1)  which  are  continuous 

from  the  left.   Then 

oo 


h(t)  =/  (v^+t-^x)^^ 


-oo 

is  a  continuous  function  of  t  which  converges  to  +cd    and  -co  as  t 
goes  to  minus  and  plus  infinity  respectively.   Hence,  there  exists  a 
value  of  t  so  that 

^x+t  =  ^x  • 

As  the  value  u  was  chosen  arbitrarily  it  follows  that  the  function 
o 

u  constructed  above  assumes  every  value  and  is  therefore  a  continuous 

X 

function  of  x.   This  completes  the  proof. 


4.  Conclusions  and  Extensions 
A  second  order  accurate  scheme  to  approximate  (l.l)  is 
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where 

Q(a,b)  =lM:^lM  +  |t.£:M±j::M.(f(a)-f(b))  . 

The  coefficients  of  the  linearized  equation  are  not  positive  no 
matter  how  small  At/Ax  is  taken.   In  general,  the  linearized  equation 
of  a  second  order  accurate  scheme  does  not  have  coefficients  of  the 
same  sign. 

Numerical  experiments  indicate  that  (4.1)  is  a  stable  approxima- 
tion to  (l.l)  and  that  traveling  waves  exist  which  propagate  with  the 
speed  s  determined  by  the  Rankine-Hugoniot  relations.   The  methods 
used  above  cannot  be  used  to  establish  the  existence  of  traveling 
waves  for  these  schemes  as  a  second  order  scheme  is  not  monotone. 

An  explicit  calculation  in  the  case  s  =  0  for  the  scheme  (4.1) 
results  in  the  proof  of  the  existence  of  multiple  traveling  waves 
which  have  limits  u  =  -1  and  u   =  +1.   In  this  case,  (2.5)  reduces 

n 

to  a  first  order  recurrence  and  as  such  defines  an  operator  /^-  deter- 
mining V  . -,  as  a  function  of  v  .   One  can  determine  explicitly  the 
x+l  X 

set  of  points  in  the  real  line  which  when  iterated  under  i\.   converge 
to  -1.   The  essentially  different  solutions  occur  from  considering 
different  branches  of  the  function  P^   . 

The  reader  should  note  that  in  general  (2.5)  defines  v^_|_j^  as  a 

function  of  v  j^w_-]>"'>^   -k+1*  ^'^   ^^^  usual  construction,  the  rela- 

2k—  1 
tion  can  be  considered  to  determine  a  function  ;^„  from  R     into 

itself.   The  point,  u  ,  in  R   ~  ,  all  of  whose  components  equal  u^  is 

a  fixed  point  of  .C  .   The  set  of  points  which  when  iterated  under  /C- 

converge  to  u  is  called  the  stable  manifold  of  ju  at  u^,  [7] • 

Lemma  2.3  is  an  explicit  parametrization  of  the  stable  manifold  of  J^ 
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at  u   for  the  class  of  monotone  difference  operators,  and  Lemma  2o4 

r>  —  1    ^ 

is  a  paraimetrization  of  the  stable  manifold  of  .'\L   at  u„«   Theorem 

1  states  that  these  two  manifolds  intersecto   With  this  approach  the 
existence  of  traveling  waves  can  be  established  when  u,  f,  and  g 
take  on  values  in  R  and  u   and  u  are  close.   The  details  will  appear 
in  a  subsequent  paper. 
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